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Two-dimensional materials have attracted increasing attention because of their particular properties and 
potential applications in next-generation nanodevices. In this work, we investigate the physical and chemical 
properties of waved graphenes/nanoribbons based on first-principles calculations. We show that waved 
graphenes are compressible up to a strain of 50% and ultra-flexible because of the vanishing in-plane 
stiffness. The conductivity of waved graphenes is reduced due to charge decoupling under high compression. 
Our analysis of pyramidalization angles predicts that the chemistry of waved graphenes can be easily 
controlled by modulating local curvatures. We further demonstrate that band gaps of armchair waved 
graphene nanoribbons decrease with the increase of compression if they are asymmetrical in geometry, 
while increase if symmetrical. For waved zigzag nanoribbons, their anti-ferromagnetic states are strongly 
enhanced by increasing compression. The versatile functions of waved graphenes enable their applications 
in multi-functional nanodevices and sensors. 

Graphene, a two dimensional monoatomic building block of a carbon allotrope, has attracted world-wide 
attention and being studied in nearly every field of science and engineering due to its exceptional charge 
transport, thermal, optical, and mechanical properties 1 " 9 and diverse potential applications 9 " 16 . In the 
electronic spectrum of pristine graphenes, valence and conduction bands cross the Fermi level at a single point in 
Brillouin zone (K space), the Dirac Point, and graphene is referred to as a zero band gap semiconductor or 
semimetal. Its particular electronic properties offer high carrier mobility and saturation velocity for future high- 
speed electronics. However, its semimetal nature also prevents its application in switching devices, such as logic 
gates and optoelectronic devices. Various methods have been introduced to open its band gap, such as functio- 
nalizing the surface, creating point/structural defects, doping other elements, and coating it on substrates 17 " 28 . 
Since ripples are difficult to avoid in two-dimensional materials and are an intrinsic feature of graphene, which 
strongly affect the electronic properties of graphene 29 " 35 , controlling the rippled structures in graphene had 
resulted in the design of functional devices 14 ' 29 , where three dimensional (3D) periodic structures were achieved 
for suspended or flexible- substrate supported graphene. The ripples in graphene can accommodate large struc- 
tural deformations, overcome adhesion failure under high levels of strains on flexible substrate, and allow the 
applications in flexible electronics. For example, the graphene with periodic ripples in one direction showed 
excellent flexible property 14 . Although theoretical studies showed that the ripples led to inhomogeneous charge 
distribution and midgap in graphene 8,36 , the effects of the ripples on the elastic, electronic, magnetic, and chemical 
properties are absent. In this work, we investigate the physical and chemical properties of graphenes and 
nanoribbons with periodic ripples in one direction, called waved graphenes/nanoribbons, based on first-princi- 
ples calculations. Our calculations show that waved graphenes are ultra-flexible upon compression. We also show 
that waved graphenes keep as semimetal, but the charge distribution is distorted. We demonstrate that their 
chemical reactivity is tunable by simply modulating the local curvatures. We further show that the band gap of 
waved armchair nanoribbon is dependent on its geometrical structure and compression, and the anti-ferromag- 
netic state of zigzag nanoribbon is enhanced by compression. 

Results and Discussion 

The waved graphene with periodic ripples along one direction, which can be experimentally achieved on a pre- 
strained substrate 14 , is constructed by creating repeated valleys and ridges along a certain direction (see Fig. 1). In 
this study, we consider two kinds of waved graphenes (WGs), one with the wave along the armchair direction 
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Figure 1 | Representative structures of (a) a flat graphene with arrows indicating directions of the periodic ripples, (b) a waved graphene with the periodic 
ripples along the armchair direction, and (c) a waved graphene with the periodic ripples along the zigzag direction. 



(ac-WG) (see Fig. lb) and another along the zigzag direction (zz- 
WG) (see Fig. lc). Waved graphenes are realized based on the calcu- 
lated lattice parameter of the graphene (C-C bond length is 1.427 A). 
The starting wavelength of waved graphenes (c 0 ) is about 25 A 
(referring to the supercell's length of the flat graphene along arm- 
chair/zigzag direction). Waved graphenes with shrunk wavelength 
(c) are realized by compression (see Fig. 2). All waved graphenes with 
various wavelengths are relaxed to investigate the elastic, electronic, 
magnetic, and chemical properties. 

Ultra-flexibility. As an indication of the stability of the waved 
graphene, its formation energy, as compared with that of the flat 
graphene (FG), is calculated as: 

E dif = (E tot (WG)-E tot (FG))/N 
where E tot (WG) is the total energy of the waved graphene under 




Figure 2 | Representative structures of flat and waved graphenes with 
various wavelengths. 



strain (s= x 100%), E tot (FG) is the total energy of the flat 

Co 

graphene (s = 0), and N is the total number of carbon atoms in 
the supercell. The calculated energy difference (E di f) between the 
waved and flat graphenes shows that the energy of the waved 
graphene linearly increases with the compression (the reduction of 
wavelength) (see Figs. 3a&d). The increment of the energy is only 
1.5% per atom with a compression up to 50% (see inset in Fig. 3a). If 
the graphene keeps flat under the compression, its energy goes 
nonlinearly up with increasing strain (see Fig. 3b and inset in 
Fig. 3d). By comparing the energy of the waved graphene with that 
of the flat one under the same strain, we see that the energy increment 
in both of them are almost equal when s < 1% (see Fig. 3c). When s 
> 1, the energy of the waved graphene is much less than that of the 
flat one (see Fig. 3), indicating that the graphene cannot keep flat 
under high compression (s > 1%). Our calculations confirm that the 
waved graphene is stable because the slope of energy increment is 
very low (—2.63 and 2.96 meV/atom for ac-WG and zz-WG, 
respectively) (see Figs. 3a&d), and can be easily obtained on a pre- 
strained substrate 14 . 

The strain energy (E s (e)), namely, the total energy at a given strain 
s minus the total energy at zero strain, is given by E di f X N. Then, we 
can obtain the tension force (F T = —dE s (s)/dc), the force constant 

(K = d 2 E s (s)/dc 2 ), and the in-plane stiffness (C = — d 2 E s (s)/d8 2 ) in 

terms of the equilibrium area of the supercell (A 0 ) 37 " 39 . The calculated 
tension forces are 0.24 and 0.41 eV/A for ac-WG and zz-WG, 
respectively. Because of the linear relationship between the strain 
energy and the compression, the force constant and in-plane stiffness 
of the waved graphene are zero, indicating that they are ultra- flexible 
and keep undamaged even the substrate is released from high pre- 
strain, which is consistent with the experimental results in literatures 
[14]. The ultra-flexibility of waved graphenes enables their applica- 
tions in stretchable devices. 

Electronic properties of waved graphenes. Similar to the lattice 
structure of waved graphenes that keeps stable under compression, 
the band structures of waved graphenes are not strongly affected. For 
the waved graphene with its ripples along the armchair direction (ac- 
WG), the calculated band structures are almost unaffected by the 
applied compression up to 50% (or the reduction of wavelength) (see 
Fig. 4). We clearly see that the Dirac cone of flat graphene is at around 
1/3 along the reciprocal direction M (0, 0.5, 0) to T (0, 0, 0) (see 
Fig. 4a), which is almost unchanged with the increasing compression 
(see Figs. 4b ~ f). We also note that more states in the conduction 
band are pushed down with the increasing compression (see Fig. 4). 
Similar to the response of ac-WGs to external strains, the graphene 
with its ripples along the zigzag direction (zz-WG) is also tolerant to 
the compression (or the reduction of wavelength). The calculated 
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Figure 3 | The calculated strain energies of ac-WG (a) and flat graphene (b) as a function of compression along the armchair direction, (c) fine-view 
comparison of the calculated strain energies of ac-WG and flat graphene with a compression up to 6%. (d) The calculated strain energies of zz-WG as a 
function of compression along the zigzag direction. Inset in (a) is the incremental energy as a function of compression. Inset in (d) is the calculated strain 
energies of flat graphene as a function of compression along the zigzag direction. 
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Figure 4 | The calculated band structures of ac-WGs under different compressions. The Fermi level is at 0 eV. 
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band structures show that the Dirac cone keeps almost unchanged 
with the increase of compression (see Fig. 5), indicating they are 
semi-metallic. The band structure also keeps unaffected under the 
compression, except the Dirac cones near the Fermi level under high 
compression, which are distorted with energy levels being pushed 
down and flattened within the reciprocal region (M to T) (see 
Figs. 5e&f). The calculated electronic structures of wave graphenes 
show that their conducting characters are unaffected by the 
compression. To further reveal the effects of compression on the 
conductivity, the charge densities of bands near the Fermi level, 
including conduction band bottom (CBB) and valence band top 
(VBT), are calculated. It is well-known that the VBT and CBB 
states of the flat graphene are attributed to the big n bonding and 
anti-bonding (tc*) states 8,9 . Similar to flat graphene, the big n and n* 
states of the zz-WG keep extended in the whole structure (see Fig. 6). 
The bending in the zz-WG occurs at the shared length between two 
hexagons (see inset in Fig. 6), which may not destroy the coupling 
between the big n bonding states of zz-WG under lower compres- 
sion. However, the coupling is distorted by high compression (see 
Fig. 5), leading to the distortion of Dirac cone. Different from zz-WG, 
the bending in the ac- WG occurs at the bridge between two hexagons 
(see inset in Fig. 7), which could result in the decoupling between the 
big n bonding states in the neighboring hexagons. As shown in 
the calculated charge densities of VBT and CBB of the ac- WG, the 
coupling between n bonding states along the wave is weakened due to 
the curvature (see Fig. 7a). The decoupling between n bonding states 
in valence band leads to the charge redistribution, which confirms 
the reduction of conductivity 8 . We see that the weak compression 
may have less effect on the conductivity of zz-WG, while the 
resistance of ac-WGs should increase. It is expected that the waved 
graphenes can be mechanical sensors and the change of resistance 
can also indicate the orientation of periodic ripples on graphene. 

Chemistry of waved graphenes. The ripples in waved graphenes 
affect not only their conducting properties, but their chemistry as 
well, due to local curvatures. According to the "jr-orbital axis vector" 
(POAV) theory 40 " 44 , carbon atoms at edges, defects, or regions of high 
local curvature are more active due to the elevated chemical potential 
and strain energy, while the carbon atoms at flat regions are relatively 
inert due to conjugation 40 " 49 . The heightened chemical reactivity of 
atoms at those sites with local conformational deformation can be 
quantitatively characterized by the POAV pyramidalization angle 40 " 44 , 
which is defined as (0 CT7r — 90°) with Q an being the common angle 
made by the Tt-orbital axis vector and the a bonds 43,50 ' 51 . The 




Figure 6 | The calculated charge densities of (a) valence band top 
and (b) conduction band bottom of zz-WG under a compression of 20%. 
The inset shows the flat graphene with arrows indicating the bending of zz- 
WG that occurs. 



calculated POAV angles of the atoms on the waved graphenes (see 
Fig. 8) show that the chemical reactivity of carbon atom depends on 
its location and the current wavelength (or compression) as well. The 
plot of POAV angles of the atoms on a waved graphene (see Fig. 2) 
appears to follow a similar wave shape as that of the graphene, which 
is different from these in a fullerene (POAV angle in C 60 is 11.64 0 ) 43 
and a single-walled carbon nanotube (for example, 5.15° in zigzag 
(10, 0) SWCNT) 49 . The Q a7Z at crest and trough deviates significantly 
from 90°, with larger POAV angles at crest (positive) and trough 
(negative), indicating the high chemical reactivity at these sites due 
to the high strain energy (see Fig. 8). The POAV angle (negative) 
increases as the location of an atom shifts from the trough to the 
middle between trough and crest, which is due to the release of 
curvature, indicating the gradual reduction of the chemical 
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Figure 7 | The calculated charge densities of (a) valence band top 
and (b) conduction band bottom of ac-WG under a compression of 20%. 
The inset shows a flat graphene with arrows indicating the bending 
of ac-WG that occurs. 

reactivity. With its location further shifting from the middle between 
the trough and the crest to the crest, the POAV angle (positive) con- 
tinue to increase and the chemical reactivity is therefore further 
enhanced. Notably, the POAV angles of atoms at crest (positive) in- 
crease and those at the trough (negative) decrease with the reduction 
of the wavelength or the increment of compression because of the 
increased local curvature (see Fig. 8), suggesting that the chemical 
reactivity can be enhanced by compressing the waved graphene. 
Considering its ultra-flexibility, the chemical reactivity of carbon 
atoms on waved graphene can be easily controlled and the chemi- 
cal reaction on waved graphene can be modulated by a simple 
mechanical compression or stretch. The waved graphenes may 
possess the similar chemistry of molecular species as those of 
SWNCT and fullerene with improved and controllable reactivity. 

Electronic and magnetic properties of waved graphene nanorib- 
bons. Different from the stable lattice and tolerant band structure of 
waved graphenes under compression, the electronic properties of 
waved nanoribbons are evidently affected by applied strains. To 



construct nanoribbons, the chirality and symmetry are considered, 
including armchair and zigzag nanoribbons with asymmetrical and 
symmetrical edges, respectively (see Supplementary Figs. SI online). 
Waved graphene nanoribbons (WG-NR) under the compression 
from 10 to 50% are directly built from the optimized structures of 
the waved graphenes under the same compression, and their edge 
dangling states are saturated by H atoms. 

For convenience, we adopt the parameter (p) to define the width of 
armchair WG-NRs, which is equal to the number of carbon dimer 
lines in the super cell (see Supplementary Figs. Sla & Sib online). 
Our calculations show that armchair WG-NRs (ac-WG-NRs) are 
direct band semiconductors and non-magnetic (see Fig. 9). The cal- 
culated band structures show that the band gaps of ac-WG-NRs with 
asymmetrical edges (p = 16) are larger than those with symmetrical 
edges (p = 17) (see Figs. 9a&g) under vanishing compression, 
because the band gap of graphene nanoribbons exhibits three distinct 
families of characteristic behaviours with width (p) and the gap 
hierarchy follows E g , 3n+1 > E g> 3n > E g> 3n+2 (where n is a positive 
integer, and p = 3n + m (m = 0, 1, or 2)) 39,52 . Interestingly, the band 
gap of asymmetrical ac-WG-NR (ac-WG-NR-a) increases at first 
with the compression up to 10% (see Figs. 9a&b), and then reduces 
with the further increment of the compression (see Figs. 9c ~ f). The 
change of band gap as a function of compression is shown in Fig. 10. 
Clearly, we see that conduction band bottom (CBB) and valence 
band top (VBT) move up and down, respectively, within the com- 
pression from 0 to 10% (see Fig. 9b), resulting in the increment of 
band gap. We also note that the energy level just above CBB (CBB-a) 
and that below VBT (VBT-b) move reversely to the movement of 
CBB and VBT under the compression. At 10% compression, the CBB 
and CBB-a, and VBT and VBT-b join together at T points, leading to 
the maximal band gap of ac-WG-NR-a. By further increasing the 
compression, CBB-a becomes new CBB and moves down, VBT-b 
becomes new VBT and moves up, resulting in the reduction of the 
band gap of ac-WG-NR-a (see Figs. 9c ~ f). Different from ac-WG- 
NR-a, the calculated band structures show that the band gap of 
symmetrical ac-WG-NR (ac-WG-NR-s) increases linearly with the 
increment of compression (see Figs. 9g ~ m & Fig. 10). We also see 
that its CBB and VBT move up and down as the compression 
increases, respectively, leading to the increment of band gap under 
compression. Although CBB-a and VBT-b also move oppositely to 
the movements of CBB and VBT, the gap spaces between CBB and 
CBB-a, and VBT and VBT-b are too large for CBB-a and VBT-b to 
overtake CBB and VBT (see Figs. 9g ~ m). The charge density of each 
energy band further confirms the movement of these bands under 
compression (see Supplementary Figs. S2 & S3 online). 
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Figure 8 | The calculated pyramidalization angles of carbon atoms within one wavelength of waved graphene with the ripples along (a) the armchair 
direction (ac-WG) and (b) the zigzag direction (zz-WG) under compression (10 ~ 50%). The insets indicate the locations of carbon atoms in 
ac-WG and zz-WG under a compression of 20%, for example. 
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Figure 9 | The calculated band structures of armchair waved graphene nanoribbons with asymmetrical edges (up row) and symmetrical edges (down 
row) under compressions. The Fermi level is at 0 eV. 



Previous theoretical studies showed that the zigzag graphene 
nanoribbons exhibit an antiferromagnetic insulating ground state 
with edge- state electrons that experience ferromagnetic coupling at 
each edge but antiferromagnetic coupling across the ribbon 52-54 . Our 
calculated exchange energies (energy difference between anti-fer- 
romagnetic and ferromagnetic states, E a f m — Ef m ) also show that 
the zigzag graphene nanoribbons without compression are anti-fer- 
romagnetic because that the energies of anti- ferromagnetic states are 
lower than those of ferromagnetic states. At the same time, we find 
that the anti- ferromagnetic ground state of asymmetrical nanorib- 
bon is more stable than that of symmetrical one because of the 
relatively higher exchange energy (see Fig. 11). Most importantly, 
we see that the anti- ferromagnetic ground states of zigzag waved 
graphene nanoribbons (zz-WG-NRs) are efficiently enhanced by 
applying compression (or the reduction of the wavelength) (see 
Fig. 11). The exchange energies of asymmetrical and symmetrical 
zz-WG-NRs (zz-WG-NR-a and zz-WG-NR-s) are increased by 
50% and 65%, respectively, under a compression up to 50%. The 
calculated band structures of zz-WG-NRs show that they are semi- 
conductors with band gaps of —0.4 eV (see Supplementary Fig. S4 
online). The calculated band structures of zigzag graphene nanorib- 
bons without compression are consistent with the literatures 52,53 , but 
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Figure 10 | The calculated band gap of armchair waved graphene 
nanoribbons with asymmetrical edges and symmetrical edges as a 
function of compression. 



their band gaps are less than those in Ref. 53 because the quasi- 
particle effect is not considered in this work. Their band gaps are 
almost unaffected by the compression, except under high compres- 
sion of 50%, which leads to the change from direct band to indirect 
band and the reduction of band gap (—15%) (see Supplementary Fig. 
S4 online). 

The dependence of compression of the electronic and magnetic 
properties of WG-NRs demonstrates that they are sensitive to 
applied strain and may be applicable to a variety of functional nano- 
devices. For example, ac-WG-NRs can be used as mechanical 
switches and sensors, and electronic devices with adjustable thresh- 
old voltage because their band gaps can be easily controlled by com- 
pression. By engineering their edge structures and applying strains, 
multi-functional nanodevices may be achievable because their band 
structures as a function of compression are also edge-dependent. 
Similarly, the magnetic response of zz-WG-NRs to applied compres- 
sion also indicates that they can be used in mechanical sensors and 
spintronics. Reversely, the edge structures of WG-NRs can be easily 
characterized by measuring their responses to compression in 
experiments. 
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Figure 1 1 | The calculated exchange energies of zigzag waved graphene 
nanoribbons with asymmetrical edges and symmetrical edges as a 
function of compression. 
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Conclusions 

In summary, we present a first-principles study on the elastic, elec- 
tronic, magnetic, and chemical properties of waved graphenes and 
nanoribbons. We find that waved graphenes show ultra-flexiblility 
and are compressible up to 50% with a starting wavelength of —25 A, 
which can be higher with a larger starting wavelength. We show that 
the conductivity of waved graphenes depends on the chirality of the 
periodic ripples. Waved graphenes with periodic ripples along the 
armchair direction show increased resistance under applied com- 
pression due to the decoupling between the n bonding states in 
valence band. The characteristics of the band structure of graphene, 
such as Dirac point, keep unaffected with the existence of periodic 
ripples, except waved graphenes with periodic ripples along zigzag 
direction under a compression of 50%, where the Dirac point is 
distorted. We predict that the chemical reactivity of waved graphene 
can be easily modulated by applying external stain to tune its wave- 
length or local curvature. The band gaps of armchair waved graphene 
nanoribbons can be controlled and the anti-ferromagnetic ground 
stated of zigzag waved graphene nanoribbons are greatly enhanced 
by the applied compression. Waved graphenes and nanoribbons may 
find applications in stretchable nanodevices, mechanical sensors and 
switches, chemical reaction, and spintronics by controlling external 
strains and engineering edge structures. Our calculations also suggest 
a simple approach to determine the edge structures of nanoribbons 
by detecting their electronic and magnetic response, and chemical 
reactivity to applied strains. 

Methods 

The first-principles calculations based on the density functional theory (DFT) 55 and 
the Perdew- Burke- Eznerhof generalized gradient approximation (PBE-GGA) 56 are 
carried out to study the electronic and elastic properties of waved graphenes. The 
projector augmented wave (PAW) scheme 57 ' 58 as incorporated in the Vienna ab initio 
simulation package (VASP) 59 is used in the study. The Monkhorst and Pack scheme of 
k point sampling is used for integration over the first Brillouin zone 60 . A 15 X 1 X 1 
grid for k-point sampling for geometry optimization and calculations of density of 
states, respectively, and an energy cut-off of 500 eV are consistently used for the 
nanoribbons in our calculations. Good convergence is obtained with these parameters 
and the total energy was converged to 2.0 X 10" 5 eV/atom. 
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